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Abstract 

The T]'g*g^ effective vertex function is calculated in the QCD hard-scattering approach, 
taking into account the 7/-meson mass. We work in the approximation in which only one 
non-leading Gegenbauer moment for both the quark- ant iquark and the gluonic light-cone 
distribution amplitudes for the ?/-meson is kept. The vertex function with one off-shell 
gluon is shown to have the form (valid for \qf \ > m^,) F rj i g » g {q\, 0, m^,) = m^,H(qf)/ (qf — 
m^,), where H(qf) is a slowly varying function, derived analytically in this paper. The 
resulting vertex function is in agreement with the phenomenologically inferred form of 
this vertex obtained from an analysis of the CLEO data on the 7/-meson energy spectrum 
in the decay T(IS') — > rj'X. We also present an interpolating formula for the vertex 
function F^ g * g (qf,0,m^,) for the space-like region of the virtuality q\, which satisfies the 
QCD anomaly normalization for on-shell gluons and the perturbative-QCD result for the 
gluon virtuality \qf\ > 2 GeV 2 . 
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1 Introduction 



We reanalyze the T]'g*g(*> vertex involving two gluons and the //-meson, in which one 
or both of the gluons can be virtual. The vfg*g^*' effective vertex function (or the rj - 
gluon transition form factor) enters in a number of decays such as J/ip — > 7/7, B — > 
(ir, p, K, K*)rf , B — > rj'Xg, T — > rj'X, T — > 7/7, and hadronic production processes, 
such as iV + N(N) — > //X, and hence is of great phenomenological importance. Its 
electromagnetic counterpart, namely the vertex 1/7*7 (equivalently the rj — 7 transition 
form factor) has been measured in the process 77* — > //, and analyzed in the perturbative 
QCD approach to exclusive processes. There is no direct experimental measurement of its 
QCD analogue, i.e., the process gg* — > rf, but inclusive decays of the heavy mesons such 
as B — > rj'X s and T(IS') — ► rj'X are probably as close as one could get phenomenologically 
to the underlying vertex. 

A first attempt to describe the rjg*g vertex in terms of the convolution of the distribu- 
tion amplitudes (DAs) of the //-meson in the lowest twist (twist-two) approximation and 
a hard-scattering kernel, involving the perturbatively calculable processes qq —* gg and 
gg — > gg, was undertaken in Ref. £Q. In an earlier paper j2j, we extended this analysis 
to the case involving two virtual gluons, i.e. the vertex rj'g*g*, and studied the effect of 
including the transverse momentum of the partons in the //-meson using the Sudakov 
formalism. Subsequent to this detailed study, the r]'g*g* vertex was studied in Refs. |3J 
and [1], ignoring the transverse momentum effects. The latter of the two included power 
corrections in 1/ Q 2n by using the running coupling constant method in the standard hard- 
scattering approach. However, both of these papers, as well as the earlier ones, neglected 
the //-meson mass effects. While ignoring the meson mass is an excellent approximation 
for the pion-photon transition form factor (equivalently the 717*7 vertex), this is not ex- 
pected to be quantitatively reliable for the case of the //-meson due to its large mass, 
rrirf = 958 MeV, in particular for the lower Q 2 region, where Q 2 is the virtuality in the 
r)'g*g(*' vertex. We present an improved calculation of the r)'g*g* effective vertex func- 
tion including the //-meson mass. In doing this, we also correct an error in our earlier 
paper due to an inappropriate choice of the projection operator of the //-meson onto 
the two-gluon state. 

A phenomenological form for the rj' — g transition form factor was proposed by Kagan 
and Petrov some time ago 0: 



, ., 2 , rn 2 n ,H(ql 0,m 2 n ,) 



HI ''V 

As an explicit form for the function H(qf, 0, m 2 ,), these authors suggested to neglect the in- 
dependence and approximate this function as a constant, with H(q 2 ,0,m 2 ,) ~ 1.8 GeV -1 
extracted from the data on the decay J /if) — > 7/7 [6J. This phenomenological transition 
form factor is in qualitative agreement with the //'-meson energy spectrum near the upper 
end of the spectrum in the process T(1S) — * rj'X measured recently by the CLEO 
collaboration 0, and also analyzed by us in a recent paper j^j. We show in this paper 
that the form of the 7/ — g transition form factor is obtained in the perturbative hard- 
scattering approach, if a light-cone wave-function for the //'-meson is used including the 
//'-meson mass. 
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The validity of the perturbative QCD formalism is expected to hold only above a 
certain gluon virtuality, which for the time-like region of q\ must be definitely larger 
than rrvii to avoid the pole. For the space-like region, there is no a pole for any value 
of q\ and also the effects of including the transverse momentum of the partons in the rf- 
meson wave-function are numerically not important, as shown in Ref. |2j. However, also 
for this case, for low values of the gluon virtuality important non-perturbative effects are 
present and the perturbative QCD formalism for the r]'g*g vertex is no longer applicable. 
In our numerical calculations, we shall set Ql = /Xq = 2 GeV 2 , where Q and /i are 
the threshold in the evaluation of the perturbative QCD kernel and the starting point of 
the evolution of the 7/-meson distribution amplitudes, respectively. For on-shell gluons, 
the rj'g*g vertex function F v > gg (Q, 0, rah) is determined by the QCD anomaly which is 
a non-perturbative result. We propose an interpolating formula for the vertex function 
F v >g* g (qf, 0, m 2 ,), which interpolates between the normalization of this function for on-shell 
gluons, as determined from the QCD anomaly, and the perturbative QCD hard-scattering 
result for space-like gluon virtualities with Q 2 > Ql where Q 2 = \qf \ +m 2 ,. For the time- 
like virtuality, an interpolating formula for the vertex function remains to be worked out, 
as one must include the transverse-momentum effects, which are large at low values of q 2 , 
in addition to the non-perturbative and the 7/-meson mass effects. 

This paper is organized as follows: In sectional we present the f/-meson projection op- 
erators onto the quark- ant iquark and the gluonic states, taking into account the 7/-meson 
mass effects, and the leading-twist distribution amplitudes for the rf- meson. In section 
we derive the rj'g*g* effective vertex function in the perturbative QCD approach and dis- 
cuss the region of its applicability. Numerical analysis of the effective vertex function is 
presented in section HJ where we have used the constraints on the Gegenbauer coefficients 
obtained by us [9 j from the analysis of the data on the T(lS') — > rj'X decay, combined 
with the corresponding constraints from the data on the process 77* — > 7/ presented in 
Ref. P]. In section we give an interpolating formula for the rj'g*g vertex function for the 
space-like region of the gluon virtuality. We conclude with a short summary in section |BJ 



2 The ^'-Meson Projection Operators 
2.1 Projection onto the Quark- Antiquark State 

The ?/-meson contains both quark-antiquark and gluonic components. Being a pseu- 
doscalar meson, its quark content can be described by the matrix element of the bilocal 
axial-vector operator in the SU(3) flavour- singlet state: 

0g(x,y) = -^=^(x) 7m 75 [x,y] ¥(y), (2) 

where ty(x) = (u(x), d(x), s(x)) is the triplet of the light quark fields in the flavour space 
and Nf = 3. The summation over the Dirac, colour and flavour indices is implicitly 
assumed in this bilocal operator. The path-ordered gauge factor, 

X 

[x,y]=Vexp\ig s [ dz^(z)t B \, (3) 
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is introduced to ensure gauge invariance of the bilocal axial- vector quark operator (J2J). In 
the gauge factor (J3J), g s is the QCD coupling constant, A^(z) (B = 1, . . . , N 2 — 1) is a 
four-potential of the gluonic field, ts are the generators of the colour SU(N C ) group with 
N c — 3 being the number of the quark colours in QCD, and the integration is performed 
over the straight line connecting the points y and x. 

Taking into account the leading- twist (twist-two), twist-three and twist-four contri- 
butions, the matrix element of the operator (0) between the vacuum and the 7/-meson 
states can be presented in the form: 



(0\O®{x,-x)\rf(p)) = ih I due«<*> 



u) + 



m 



* x 2 . 



[u 



+ 



u 



mz 



u and 1 - 



2(px) 

(4) 

u are the momentum 
2u -1. The 
contain 



where p M is the four-momentum of the ?/-meson (p 2 

fractions carried by the quark and the antiquark inside the r/-meson, and £ 
function (f>!S\u) is the twist-two distribution amplitude (DA), and A(u) and M(u 7 
contributions from operators of twist-two, twist-three and twist-four. This set of the 
quark- antiquark two-particle DAs is completely analogous to the one defined in Ref. JUj 
for the case of the SU(3) flavour-octet pseudoscalar mesons. The decay constant f v i is 
defined in the local limit of the matrix element as 



(0|Og?(0,0)|T/'(p)> = (0|-^*(0) 7m 75*(0W(p)) = ifn'P,, 
Thus, the DAs satisfy the following normalization conditions: 
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(5) 



du0{u) 



duM{u) = 0. 



(6) 



The normalization condition for A(u) is not fixed by Eq. (J3J). This DA is related to the 
two- and three-particle DAs of lower twists with the help of the equations of motion (see 
Eqs. (3.2) and (6.11) of Ref. [10] ) and, thus, its form (and also the normalization) is 
implicitly fixed. 

We now proceed to work out the DAs of the 7/-meson on the light-cone. To that end, 
we construct two light-like vectors and P M (z 2 = and P 2 = 0) of the same dimensions 
as the two vectors x M and p^ at hand: 



x 



A 1 2 



- (px) — (px) 2 — m 2 ,x 2 



(7) 



(px) + yj(px) 2 — m 2 ,x 2 



Their scalar product has a rather complicated form: 

2 [(px) 2 



(Pz) 



m 2 ,x 2 ] 



(px) + J (px) 2 — m^iX 2 



(9) 
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which, however, can be simplified if one assumes that the separation between the quark 



and the antiquark inside the 7/-meson is light-like (x^ 
light-like momentum P^ introduced in Eq. (JHJ) then takes a simple form: 



and x z 



2(pz) ' 



(Pz) = (pz) = (px). 



0). The 



(10) 



In terms of such light-like vectors, one can give definitions of the DAs of the 7/-meson on 
the light-cone (similar to the ones for the 7r-meson [10J) following from the matrix element 
of the axial- vector current: 



(Q\0<$(z,-z)W(p)) = if„ I due^ 



P„ 



u) + 



m v fZ » (g) 
2(Pz) yr >' 



u 



;n) 



Comparing it with Eq. (j3J), one obtains the relation: 

B(u) = <#>(«)- <* (ff) 



u 



implying that the Lorentz invariant amplitudes <fity(u) and M(u) can be interpreted as 
the ?/-meson DAs of the nonlocal axial-vector operator at a strictly light-like separation. 
It is easy to see the usual normalization conditions of the light-cone DAs: 



(12) 



dug<*\u) = l. 



(13) 



Two more bilocal matrix elements define DAs of the twist-three: 



(o|- 



1 

:*(x) 75 [a;,-x]*(-x)|V(p)) = -ifrfH j due*W<j>«\u), (14) 



(o|- 



-.^{x)a a pj 5 [x, -x]$(-x)\ri'(p)) 



(15) 



x {p a x p - pfsx a ) / du e^^iffl) {u 



where = m^,/(2m 11 + 2m ( j + 2m s ) is a chirally enhanced factor including the masses m q 
(q = u, d, s) of the light quarks. Note that after the transition to the light-cone vectors P M 
and z^, the forms of Eqs. ([14)1 and (|15J) are not changed in the twist-three approximation 
and, thus, the Lorentz invariant twist-three DAs coincide with the light-cone twist-three 
DAs. There also exist twist-three DAs defined by the matrix elements of the three-particle 
(quark-antiquark-gluon) operators. For the light pseudoscalar flavour-octet mesons, it has 
been shown in Ref. JTU] that their contributions are numerically small. We expect that 
this also holds for the SU(3) flavour-singlet meson and neglect such contributions. In 
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general, one should also take into account the twist-four DAs from the matrix elements of 
three-particle nonlocal operators. The self-consistent set of such DAs was introduced in 
Ref. |11[IT2] for light pseudoscalar flavour-octet mesons; their explicit updated forms can 
be found in Ref. ^Oj- While the effect of the chirally enhanced bilocal matrix elements 
of higher twist operators is an interesting issue to study, we neglect here the twist-three 
and twist-four contributions to the //-meson light-cone wave-function. 

Starting from Eq. (fTTj) . one can calculate the light-cone projection operator of the 
7/-meson onto the quark- ant iquark state in the leading-twist (twist-two) approximation 3 : 

1 

(0\*U*)*Jf*(-*)W<P)) = ^ [TtiPlHii&P* -7= / due*M$\u), (16) 







where the two sets of indices (iaa) and (J (3b) describe the components of the quark ^f(z) 
and antiquark z) fields in the Dirac, colour and flavour spaces, respectively, and 
(P7) = -P M 7^. After the Fourier transformation of Eq. (|16|) to the momentum space (see 
Eqs. (A. 4) and (A. 5) in Ref. [3]), the result can be presented in the form: 

//V [75^7)]^* 4=$ } («) = 2(np) / ^e^H0\*U*)*M-*)W(p))> ( 17 ) 



where z^ = z~n^, and is a dimensionless light-like vector (n 2 = 0). In the momentum 
space, the quantity 



can be interpreted as the ?/-meson light-cone projection operator onto the state of the 
incoming quark and antiquark The quark- ant iquark twist-two light-cone DA, (f)^,(u), 
can be defined as follows: 

ifrf$\u) = 2 J ^-e-^{0\n^(z,-z)\rj f (p)) , (19) 

where the bilocal axial- vector operator is defined in Eq. Some words about the 

vector n M are in order. While well-defined in the position space through the light-like 
separation between the constituents inside the ?/-meson, this vector is arbitrary in the 
momentum space and its explicit choice is dictated by the process at hand. The most 
natural choice of this vector is the four-momentum of some massless particle which appears 
in a given process (see, for example, Ref. [SUED- If such a vector is absent, the light-like 
vector n M should be constructed from existing four-vectors explicitly. A representation of 
the vector will be given in the next section. 

3 The path-ordered gauge factor J3J is assumed to be between the quark fields and is suppressed to 
simplify the presentation. 
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2.2 Projection onto the Gluonic State 



The gluonic content of the ?/-meson can be described with the help of the partially 
traceless and symmetric bilocal gluonic operator 4 : 



G, a (x)[x,y]G u a (y) +G va (x)[x,y]Gf(y)] - 9 -f G^(x)[x,y]G a %y), 



(20) 

where G^ u (x) and G^ u (x) are the gluonic field strength tensor and its dual tensor, re- 



spectively, and the path-ordered gauge factor [x, y] should be taken here in the adjoint 
representation of the colour SU(N C ) group. In the above operator, the summation over 
the colour indices is implicitly assumed. Starting from the gluonic operator ([20)1 . one can 
introduce the gluonic twist-two, twist-three, and twist-four DAs, similar to the quark- 
antiquark DAs corresponding to the axial- vector bilocal operator (J2J), as follows: 



(0|<5jg(a;, -x)\r]'(p)) 



frf Gj, 



du e^ px) 



+ 



2JN f 
v o 

(px) 



nr 



P^Pv 



(pi g ,\u) + m 2 ,x 2 C{u) 



nr. 



4 



B(u) 



(21) 



where C F = (N 2 - 1)/(2N C ). 

In the local limit, the bilocal gluonic operator (|2*0)l vanishes [Oj$(0, 0) = 0] due to the 
following operator relation: 



G^(0)G/(0) = ^G^(0)G^(0), 



'.22) 



which gives the normalization conditions for the Lorentz-invariant gluonic DAs introduced 
in Eq. (jUJ): 



dud$(u) 



0. 



duB(u) = 0. 



(23) 



These conditions leave an arbitrariness in the choice of a constant prefactor in the matrix 
element f!21|) . The motivation of the choice made in this paper will be commented in the 
next subsection. The normalization of the DA C(u) is not determined in the local limit 
of the operator (|20|1. but this DA can be related with the other two- and three-particle 
DAs of lower twists using the equations of motion in the same manner as has been done 
for the flavour-octet pseudoscalar mesons [T01HIJH2] " 

Following the discussion in the preceding subsection, let us assume that the two gluons 



are separated by the light-like vector x. 



0). Thus, the matrix element of 



the bilocal gluonic operator can be presented in terms of the gluonic light-cone DAs: 

i 



2JN 



u 



(24) 



o 



4 The method of construction of the completely traceless quark-antiquark and gluonic operators can 
be found in Refs. JH] and 
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where the four- vector P M is defined in Eq. (JTUJ). Comparison of this equation with Eq. (|2~T|) 
leads to the following relation: 

B(u)=g%\u)-$\u). (25) 

This shows that the Lorentz-invariant DAs, <ply\u) and O(w), are connected with the 
light-cone DAs of the //-meson defined by the bilocal operator (j2*U|) with a strictly light- 
like separation between the gluon fields, in the same manner as for the quark-antiquark 
DAs, 4>ff(u) and B(-u), defined by the bilocal axial- vector operator @. The normalization 
conditions for the gluonic light-cone DAs are: 

•Mf„,\ _ n I m, Ja) 



du^\u) = 0, J dug^\u) = 0. (26) 
o 

As we restricted ourselves to the leading-twist part of the quark-antiquark component of 
the 7/-meson only, we also do not consider any further the higher-twist DAs in the gluonic 
component of the ?/-meson. 

The gauge-invariant definition of the gluonic light-cone DAs (124}) is given in terms of 
the gluonic field strength tensor G^ v {z) and its dual G^ u (—z), though the usual Feynman 
rules involve the gluonic four-potential A^(z). A possible way to get the required matrix 
element is to use the relation between the field strength tensor and the four-potential in 
the light-cone gauge (n^A^x; n) = 0) [T7j : 

oo 

A*(x; n) = n" j G%{x + an) da, (27) 
o 

so that it can be applied to Eq. (}2"4"j) after its contraction with the light-like vectors n M 
and n v . In the leading-twist approximation, the projection operator of the //-meson onto 
the gluonic state in terms of the gluonic four-potentials can be written as follows [3] 5 : 

<0|^Ki(-.)|V«> ^^-^f ) ^>|^, (28) 

V Q 

where AHx)A^(y) = [A^(x)A^ (y) — A„(x)A^(y)]/2 is the bilocal operator antisym- 
metrized in the Lorentz indices. The gluonic matrix element (I28j) contains the Lorentz 
structure e^pa z p P a /(zP). To perform the Fourier transform to the momentum space, 
it is convenient to introduce the dimensionless light-like vector n a (its explicit form in a 
specific frame is given in Eq. (A. 3) of Ref. |Hj) so that z a = z~n a . The light-like vector P^ 
defined by Eq. becomes independent of the variable z~: 



'Starting from here we suppress the path-ordered gauge factor [z, —z] in the matrix elements. 



and the Lorentz structure under consideration transforms as 



- pupa 



z ppa 



- pu per 



n ppa 



' pupa 



n p p a 
(np) 



(30) 



Note the invariance of this Lorentz structure under the transformations P„ — > P' 



Pu + Cn u or n b 



n, 



n^ + C Pfj,, where C and C are arbitrary factors. The last 



expression in Eq. (J3Uj) corresponds to the specific choice of the first factor: C = m 2 ,/2(np). 
The invariance of the Lorentz structure (|30J) under the transformation of the light-like 
vector was pointed out in Ref. [3] for the case of the massless n'- meson (i.e., with 
p 2 = and P M = 

After performing the Fourier transformation of Eq. (|28|) . we obtain the following result: 



frf&AB £pup*nPp° $\u) 



8N C ^/N} (np) u{l-u) 
In the momentum space, the quantity 

>(9) _ 



2(np) 



dz 
~2^ 



^<0|A£(z)^(-z)|i/(p)>. (31) 



V, 



&pvp<jnPp 



pA;uB 



AN C ^/W f (np) 



(32) 



is the projection operator of the r/-meson onto the state of two incoming gluons ^3]. The 
gluonic twist-two light-cone DA can be defined through the following matrix element: 



U^(u) 



f 



dz 



*M(0\ n >>G tux {z)n v &' a {-z)\r/{p)), 



(33) 



Cp(np) J 2n 

where the summation over the colour indices is implicitly assumed on the right-hand side. 

2.3 The ?/-Meson Light-Cone Distribution Amplitudes 

The leading-twist light-cone DAs of the 7/-meson can be presented as infinite series in the 
Gegenbauer polynomials Cn 2 (u — u) for the quark-antiquark component and Clj\(u — u) 
for the gluonic one [1311211201121111211221: 



1+ £ A n (Q 2 )Cn 



u — u) 



even n>2 



<f>^(u,Q 2 ) = u 2 u 2 B n (Q 2 ) C 5 J_\(u - u) 



5/2 



(34) 
(35) 



even n>2 

where u — 1 — u, and the following notations are introduced for the Gegenbauer moments: 



A n (Q 2 ) 
B n (Q 2 ) 



a s ([i, 



7+ 



Pn D n KPO 



a s (Q 2 



a s (Q 2 )\ 
8 



+ P ( n 9) B^(pl) 



a s (Q 2 ) 

l^(Q 2 )\ 



(36) 
(37) 



These equations contain the quantities 7™, defined as: 

1 



1QQ + iGG ± J {T QQ - 1 GG ) 2 + ^QGlGQ 



(3* 



where the anomalous dimensions 7™.- = Q,G) are 



7qq 



7q G 



7gq 



7^g 



A) 



3 + 



(n+l)(n + 2) 



n+l 



n(n + 3) 



/3b 3(n + l)(n + 2) 
Po (n + l)(n + 2)' 



(39) 



(n + l)(n + 2) 



n+l ' 



+ 1. 



Here, /3o = (HiV c — 2n/)/3 is the first coefficient in the expansion of the QCD /5-function 
and rif is the number of quarks with masses less than the energy scale Q entering in the 
DAs ()34|) and (|33j) . Note the difference between the quark numbers Nf = 3, used earlier 
in the context of the wave-function of the r/-meson, and n/ in the /3-function; the former 
reflects the quark content of the SU(3) flavour-singlet meson, while the latter is the number 
of the active quarks, which (together with the gluons) determine the renormalization 
effects in the LCD As of the rj'- meson, probed at a virtuality Q 2 . The discussion of these 
anomalous dimensions in the one- and two- loop accuracy can be found in Ref. [22] ■ Our 
choice of the constant prefactor in the gluonic matrix element (j21J) is reflected in the 
non-diagonal anomalous dimensions presented above. 

The anomalous dimensions (|3*^|) also define the quantities and pffl entering in the 
Gegenbauer moments and (j3T|) as follows: 



p (9) 
In 



6 



7+ Iqq 
^QG 



1 n 



Tqg 



6 7™ 



7qq 



(40) 



The prefactor in the gluonic light-cone DA (|3*3|) is chosen as u 2 u 2 according to Ref. 
which accounts for the appearance of the factor uu in the denominator of the integral 
function in the gluonic matrix element ()28|) . Such a functional dependence is also in 
agreement with the general expression for the "asymptotic" form of the LCDAs resulting 
from the conformal symmetry [TiH ll2|. 

We use an approximate form for the ?/-meson light-cone wave-function in which only 
the first non- asymptotic term in both the quark- ant iquark and gluonic DAs is kept. Thus, 



u, Q 2 ) = 6uu [l + 6(1 - 5uu) A 2 {Q 2 )] 



u,Q 2 



5u 2 u 2 (u-u)B 2 (Q 2 ), (41) 



where the explicit forms of the Gegenbauer polynomials C, 



3/2 



u 



u 



and C l 



5/2 



u — u) 



have been used. The free parameters B^\/j,q) and B^\ixq) (the Gegenbauer coefficients), 
which enter in the Gegenbauer moments ^(Q 2 ) and ^(Q 2 ), are not determined from 
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first principles and have to be modeled or extracted from a phenomenological analysis of 
experimental data. We defer the quantitative discussion of these Gegenbauer coefficients 
to Sec. H 

There is an additional point concerning the scale The usual choice of the scale 
— 1 GeV 2 made in the analysis of the pion form factor and the 7r° — 7 transition form 
factor represents the initial value for the evolution of the LCDAs. As the mass of the rf- 
meson is of order 1 GeV, a more realistic choice of the parameter lIq is /i 2 , = 2 GeV 2 , which 
we shall adopt. Consistent with this assumption, and with the choice Q 2 = \q 2 \ + m 2 ,, 
where q 2 is the total gluon virtuality in the r]'g*g* vertex, we shall also set Ql = 2 GeV 2 
in the calculation of the perturbative kernel. 

3 The rj'g*g* Effective Vertex Function 

In the momentum space, the effective rj'g*g* vertex can be extracted from the invariant 
matrix element of the process 7/ — > g*g* with the help of the relation 6 : 

M = MM+M® = F v/9 ^(q 2 ,qlm 2 vl )5 AB E^ef;Eg*q lp q 2(7 , (42) 

where Ai^ and Ai^ 9 ' are the contributions from the quark- ant iquark and gluonic compo- 
nents of the //-meson, respectively, and and ef (i — 1, 2) are the four- momenta and the 
polarization vectors of the final virtual gluons. The four-momentum of the //-meson is re- 
lated to the four-momenta of the gluons by energy-momentum conservation: p p = gi^+q^- 
The individual contributions Ai^ and Ai^ to the invariant amplitude (}4*2*j) can be cal- 
culated by using the //-meson projection operators onto the quark-antiquark (jl8j) and the 
two-gluonic (J32*|) states, yielding 

1 

= ifo j du <$\u, Q 2 ) V$ b . iaa S ab [T^f, (43) 


M " = f /^^^X^fe, (44) 



where and are the quark-antiquark and gluonic hard-scattering kernels calculated 
in the perturbative QCD, respectively. The factor 1/2 in Eq. f)44j) takes into account the 
two identical gluons in the //-meson. To go further, it is necessary to define the light-like 
vector rifj. appearing in the individual invariant amplitudes ()43|) and ()44j) in terms of the 
physical vectors of the problem under study. 

3.1 Specifying the Light-Like Vector n M 

We now proceed to express the vector n p in terms of the gluon momenta and q 2p in 
the T]'g*g* vertex. In particular, we consider the case when both the gluons are off the 

6 The difference in the phase factor i between the definition of the effective T)'g*g* vertex function in 
this paper and in Ref. is related to the corresponding difference in the definitions of the projection 
operators of the ry'-meson onto the quark-antiquark and gluonic states. 
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mass shell and have virtualities q\ and q\ comparable to the r/-meson mass squared. 

Let us assume, to be definite, that the virtualities of both gluons are time-like (qf > 
and q\ > 0). In that case, the two light-like vectors nj* and can be constructed 
from q lfl and q 2fl (or, equivalently, from = q lfl + q 2/I and q 2fl ): 



n< t ] = C n ^q 2 Pn + qi-qi- ™ 2 , ± A K (m v , ^Jq 2 , y^gf) g 2M | 



(45) 



where C n is an arbitrary factor and Ak(o, b, c) is a kinematic function, defined as follows: 
Ak(o, b, c) = \/ (a + b + c) (a — b — c) (a — b + c) (a + 6 — c). (46) 

In the limit of neglecting the ?/-meson mass (i.e., setting = 0), these vectors are 
simplified to: 

71^ = C n {2q\ V , + \q\ - q\ ± \q\ - q 2 2 \] q 2ll ) , (47) 
as Ak(0, ygf, a/*?I) = l<Zi ~~ Thus, if q\ > gf, then these vectors have the forms: 

- 2C n (g 2 2 g^ + g?<fc M ) , (48) 
2C„ g 2 p M , 



(-) 



while for g 2 < gf they are: 



(-) 
v- 

(+) 



If we take = n]L (n M 

,2 ^ „2 /„2 



= 2C n + q\q 2i ) , 

= 2C n q 2 p^. 

then in the limit m v ' = 0, the light-like vector P M 



(49) 



vanishes for q[ < q\ (qf > q 2 ) and, hence, the projection operators onto the quark- 
antiquark (JTHJ) and gluonic ((H2} states also vanish. To get non-vanishing projection oper- 
ators in both regions, the vector should be taken as: 



n» = n^<d{ql-ql)+n^<d{ql-ql), 



(50) 



where Q(x) is the unit step function. To write this vector uniformly, it is convenient 
to introduce the total gluon virtuality q 2 , the asymmetry parameter uj, and the relative 
//-meson mass squared rj as follows: 



2 2,2 

q = ?i + 9 2 > 



(51) 



If we also redefine the factor C n as C n = C n /q 2 , the vector (|5Uj) can be rewritten in 
the form: 



C n [(1 - + (w - ?/ + wA) g 2/ J 



where the function A can be written as 



7. 



A=,/l-H? + ^. 



(52) 



(53) 



7 Note, Ak and A are different functions. 
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In the massless limit of the rf- meson (rj = 0), A = 1. With this specification of the 
light-like vector n^, its scalar product with the ?/-meson four-momentum p^ is: 

q 2 ~ 

(pn) = -— C n u\(uj - r] + u\) , (54) 
which allows us to write the light-like vector P^ (|29|) in the following form: 

Pm = 2uj\ i( u ~V + UJ ^)P^ + 2 VQ2 t ,]- (55) 



With this, the Lorentz structure in the projection operator (J32J) is reduced to the form: 

(56) 



(np) ujX q 2 A q\ — q\ 



Note that both the vector P^, defining the Dirac structure of the rj'-meson projection 
operator onto the quark-antiquark state (j!8j) . and the Lorentz structure considered above, 
which comes from the projection operator onto the gluonic state (j!8|) . are independent of 
the choice of the factor C n . 

For the case when both the gluons have space-like virtualities (qf < and q\ < 0), 
the same set of quantities can be introduced as in Eq. (JoTj) . with the obvious difference 
that the total gluon virtuality and the relative ?/-meson mass squared are negative. This 
difference does not change the result obtained for the time-like gluon virtualities and, 
thus, Eqs. (J52J) . ()55|) and (j56|) are valid in this case also. 



3.2 The Quark Part of the r]'g*g* Vertex 

The light-like vector P^, and hence the ?/-meson projection operator ()18j) onto the quark- 
antiquark state, is completely defined in terms of the physical vectors - the four-momenta 
of the gluons. Hence, we can calculate the quark part of the r]'g*g* effective vertex function 
starting from the invariant amplitude ()43|) . obtaining the hard-scattering kernel: 

[T ( 3 ) r /3 = y*P;AB {u ^ - p; _ qu _ q2) £ M e B^ (57) 

where the Dirac indices i and j are not shown explicitly; the exact expression for the 
effective quark-antiquark-gluon-gluon vertex in the lowest order in perturbative QCD can 
be found in Eq. (3.3) of Ref. [2]. This yields the following result: 



where Q 2 = \q 2 \+ m 2 , and the summation over the quark colour and flavour indices have 
been performed. This matrix element contains the contributions from both diagrams 
presented in Fig. Q as the contribution from the second diagram can be transformed to 
the form of the first one with the help of the symmetry property of the quark-antiquark 
light-cone DA: 4>^(u, Q 2 ) = <fi^(v,, Q 2 ). It is easy to see that the parameter C defined by 
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Figure 1: Leading Feynman diagrams contributing to the quark part of the rj'g*g* vertex. 



Eq. (2.3) of Ref. [2J is connected with the 7/-meson decay constant f v i used here by the 
following relation: C = ^/Nj f v ' ■ 

Taking all this into account as well as the Ansatz (|42|) for extracting the T]'g*g* effective 
vertex function, the result for the quark part can be written as follows: 



4vra s (Q 2 ) Uy/N} 



q 2 N c uX 



du^(u,Q 2 



uj{\ + A) + rj(u — u) 
1 + uj{u — u) — 2uur\ + ie 



(59) 



The quark-antiquark contribution, in the approximation of keeping the first two terms in 
the corresponding distribution amplitude 4>^(u,Q 2 ) [see Eq. f)41|l ]. can be conveniently 
written in the following form (similar to Eq. (3.7) of Ref. j2j): 



m 2 , A 



G^\u,r ] ) + 6A 2 (Q 2 )G ( 2 q \u;,r ] )), (60) 



where A 2 (Q 2 ) is the Gegenbauer moment defined by Eq. (J36|) . and the functions G^\u, 77) 
and G 2 \oo, if) are: 



G^(u, V ) 



10 



1 - A + 

2r/ \ 77 

1-^(1 -A) 

77 



l--(l-A) 
V 



-Lin 

2uo 



l+LU 



l-OJ 

(l-A) 2 -fl-^) (1-A) 



+ A J(u,rj) 



15/ 

5 77 



(61) 
(62) 



2uj 



l + LU 



l-UJ 



+ A J(u;,?y) 



The results given above are presented using the function: 

1 



du 



1 + to{u — u) — 2uw] + ie ' 



(63) 



whose explicit form and the asymptotic behaviour can be found in Appendix B of Ref. [2] . 
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Figure 2: The functions fo(u), f 2(^)1 and ^(^) describing the large \q 2 \ asymptotic 
behaviour of the r]'g*g* effective vertex function, with the virtualities of the gluons having 
same signs. 



When the ?/-meson mass is negligible in comparison with the total gluon virtuality \q 2 \ 
(i.e., in the limit of small 77), Eq. (JoTij) reduces to: 



^a s (\q 2 \) 3f„y/N f 



N r . 



{/oM + 6A 2 (|g 2 |)/ 2 M} 



(64) 



(in accordance with Eq. (3.10) of Ref. [2]) as the functions Gq(u,tj) and G 2 g \uJ,r]) pre- 
sented above are dominated by the term G[ q \u,r]) — rj fi(u) in this limit, and the func- 
tions fi(uj) (i = 0, 2) are defined as follows: 



/oM = — 
/2H 



In 



2u I 

,2\ 



10 



12u< 



[3(5-cu 2 )/ M-10]. 



(65) 
(66) 



The dependence of these functions on the asymmetry parameter u is presented in Fig. |2j 
The asymmetry parameter is defined in the interval — 1 < uj < 1, but as the func- 
tions /o(u;) and f2{u) are symmetric, it is sufficient to present them in the region of 
positive values of the argument. 

If one of the gluons is on the mass shell, for example, the second one (g| = 0), the 
asymmetry parameter u; = 1, A = 1 — 77, and it is easy to get from Eq. ()59|) the following 
result: 



4ira s (Q 2 ) 



1 1 
q{ - m z , 



(67) 
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Figure 3: Leading order contribution to the gluonic part of the rj'g*g* vertex. 



In the approximation of the r/-meson DAs adopted here, the first inverse moment of the 
quark- ant iquark twist-two DA is: 

1 

/ ^y(H,Q 2 )=3[l + i 2 (Q 2 )]. (68) 
o 

Thus, in this case, the rj'g*g effective vertex function is defined by the first inverse moment 
of the f/-meson quark-antiquark DA and it has a pole form, with the pole being at 
q\ = m 2 ,. Of course, this result is not supposed to be used in the vicinity of the singular 
point q\ ~ m 2 ,; the perturbative form sets in at a higher value of qf, typically q\ ~ 2 GeV 2 . 
Also, as discussed in Ref. [2], this pole behaviour results from ignoring the transverse 
momentum in the definition of the //-meson wave-function and, hence, is not physical. 



3.3 The Gluonic Part of the rfg*g* Vertex 

The gluonic part of the rfg*g* vertex function originates from the diagrams presented in 
Fig. El To calculate this effective vertex function with the help of the invariant ampli- 
tude (|44jl . the corresponding hard-scattering kernel T^f* can be obtained from the effective 
four-gluon vertex V^^ D (qi, ft, ft, ft), which is defined in Eq. (4.3) of Ref. [2], as follows: 

[Tk% D = V^ D (- gi , -ft, up, up) efcg*. (69) 

Substituting Eq. (J51)Jl into Eq. (J3~2"|) . the //-meson projection operator onto the two-gluon 
state, the invariant matrix element (1441) can be rewritten in the from: 

It - <?! 



= !*=^ [ du ^ Kg2) £P t 9lA ? T 8cd [T^C, (70) 

AN cX /m\ J uu q\ -q\ K ' 

v 



where the parameter A is defined in Eq. (J53|) . Comparison of this matrix element with 
the one given in Eq. (4.1) of Ref. [2], with C = f V '\/N]f, shows that the two expressions 
differ by the factor Q 2 /[2N f X(ql - ql)} = Q 2 /(2N f Xujq 2 ). Note also the difference in the 
factor i in the definitions of the gluonic projection operators in this paper and in Ref. |2j 
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(this difference then also reflects itself in the definition of the quark- ant iquark projection 
operator). We have now understood this mismatch, related to two errors made in Ref. [2]: 
First, the factor 1/2 is due to the identity of gluons in the r/-meson, which was missed 
in Ref. [2]; Second, the factor Q 2 /(q 2 — ?|) is required by the Bose symmetry of the final 
gluons in the process rf — > g*g* described by this amplitude, also overlooked in Ref. [2]. 
Finally, the parameter A in the denominator of A4^ enters as we now take into account 
the ?/-meson mass; A = 1 for the case of the massless ?/-meson. Taking into account 
the difference between Eq. ([TO]) in this paper and Eq. (4.1) of Ref. [2| pointed above, the 
corrected result for the gluonic part of the rj'g*g* effective vertex function can be obtained 
from Eqs. (4.7) and (4.8) of Ref. |2j, which now reads as follows: 



4na s (Q 2 ) 5/, 



'/' 



B 2 {Q 2 ) G ( i\uj,r]) 



i)'g*g* 



2^/N 



[71) 



f 



where the Gegenbauer moment .82 (Q 2 ) is defined in Eq. ()41j) . and the function G < 2 s \uJ,r]) 
has the form: 



G { 2 9 \uj, V ) 



— / auuu(u 



UJ 



u 



rj + uj (u — u) 



1 + UJ (u 



u 







2uur] + ie 



(72) 



5 2 
^ + - 

6 7] 



4uj 2 



+ r] 



2 
V 



f + 2uj 
2 + uj 2 







" 4uj 2 ' 


In 


1+UJ 


V-- 






v . 




rf 




1-UJ 



+ 



iuj' 



rf 



Auj 4 
~rf~ 



J(UJ,T)). 



This function is symmetric in its first argument under the change uj — > —uj: G| 9 (— OJ, rj) = 
G^\uj, rj), in accordance with the requirement of the Bose symmetry for the rj'g*g* vertex. 

In the limit of the large total gluon virtuality (|g 2 | 3> m 2 ,), the gluonic part of the 
r]'g*g* effective vertex function simplifies and can be expressed as follows: 



F (9) 
r)'g*g 



(q 2 ,uj,0) 



47ra s (|g 2 |) 5f v , 



2,/N 



B 2 (\q 2 \) 



where the function g2{uj) has the form: 



3/o(o;) 



6uj 2 



(73) 



(74) 



Here, fo(uj) is the function defined in Eq. (|B3)l. Note that the function g2{uj) is symmetric 
under the exchange uj — > —uj: g2(—uj) = g2{u), in agreement with the observation made 
in Ref. j^j. At uj = it has a non- vanishing value: #2(0) = 1/15, and at the end 
points, uj = ±1, it has the value g 2 (±l) = 1/6. The explicit dependence of gi{uS) on the 
asymmetry parameter uj is presented in Fig. El 

If one of the gluons is on the mass shell, say, the second gluon (g| = 0), the gluonic 
part of the r]'g*g effective vertex function can be written in the following form: 



47ra s (Q 2 ) 5/„, 



2 2 
q{ - tntj 



2,/N 



B2(Q 2 )G i 2 9 \l,r ] ), 



(75) 
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Figure 4: The real (solid curve) and imaginary (long-dashed curve) parts of the func- 
tion G?2 : (1,77), where 77 = rrvL/qf, as a function of the momentum squared q\ of the 
virtual gluon. The vertical short-dashed lines, drawn at q\ = — 1 GeV 2 and q\ — 2 GeV 2 , 
demarcate the regions of applicability of the perturbative results, which are to the left 
and right of these lines for the space-like and time-like gluon virtualities, respectively. 



where r\ = m 2 ,/qf, the value A 



1 — 77 was taken into account at uo = 1, and the function 



G { 2 9 \l,v) is: 



G 



(9), 



5 

•3?/ 



2 


4 


1 


I 




4 " 








1 - - 






rf 


77 3 


V 


V. 




7] z 



ln(l — 77). 



(76) 



The dependence of this function on the gluon virtuality q\ is presented in Fig. with 
the value G^(1,0) = 1/6 corresponding to the large q\ asymptotics. In line with the 
quark- ant iquark part (J67)) of the rj'g*g effective vertex function, the gluonic part (J75j) has 
the same pole behaviour at q\ = m 2 ,. 

To conclude this section, we present the expression for the rj'g*g* effective vertex func- 
tion resulting from the perturbative QCD analysis. Combining the quark- ant iquark (|60p 
and the gluonic (fTTj) parts of the rj'g*g* effective vertex function, the perturbative result 
for the vertex can be expressed as follows: 

4vra s (Q 2 ) 



F 1 * * 

1 v /g*g* 



2 2 2 



v3/V 



(77) 



uj,r}) 



6A 2 (Q 2 )G ( 2 q) (uj, v )-^B 2 (Q 2 )G 2 
6 



(9), 
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where the total gluon virtuality is q 2 
side are defined in Eq. (|5T|). and Q 2 = 



-■ qf + g§, the variables and 77 on the right-hand 
\q 2 \ +rn 2 ,. Being a perturbative-QCD result, this 



expression is valid in the large \q \ region. 

As mentioned above, both the quark- ant iquark (|H7jl and the gluonic (J7H|) parts of the 
i]' — g transition form factor have the phenomenological form (0) which allows to extract 
the slowly varying function H(qf, 0, m 2 ,): 



H(ql,0,m 2 ) 



±na s (Q 2 



m 



V3f v , 



l + A 2 (Q 2 



B 2 (Q*)G$\l, V ) 



(78 



The non-trivial dependence on the f/-meson mass is coming through the function G 2 f (1, rj) 
with rj = m 2 ,/q 2 , which is to be traced back to the gluonic component of the //-meson. 



4 Numerical Analysis 

As demonstrated earlier, the correction due to the mass of the f/-meson appears already 
in the leading-twist (twist-two) light-cone approximation. It is of interest to know nu- 
merically the effect of including the //-meson mass in the rj'g*g* effective vertex function. 
To work this out, we need to specify the input values for the various parameters. To that 
end, we note that the 7/-meson decay constant can be related with the flavour-singlet 
decay constant f\ ~ 1.17 f^, where f n = 133 MeV is the 7r-meson decay constant, as: 
frf — /iCos#i with the mixing angle 9\ ~ —9.2° [23] • With this, it is easy to check that 
f v > ~ 2/ 7r /y / 3- In estimations of the effective vertex function, the strong coupling a s (Q 2 ) 
is used in the two-loop approximation with the QCD scale parameter A™. = 305 MeV 
corresponding to four active flavours (rif = 4). The central values of the c- and 6-quark 
MS masses, rh c — 1.3 GeV and fh^ = 4.3 GeV [23; were used for the separation of regions 
with different active-quark flavours in the strong coupling a s (Q 2 ). In this context we also 
specify the constrained parameters ()38|) and (|4T)|) in the r/-meson DAs: 7+ ~ —0.645, 
7^ ~ —1.421, — 2.863, and p% — —0.010, calculated for nf = 4. For the starting 
point of the evolution scale we take //q = Ql = 2 GeV 2 . We do not include errors on these 
parameters, as they are relatively small. 

The largest uncertainty in the f]'g*g* effective vertex function is due to the gluonic 
Gegenbauer moment B2(Q 2 ) (j37|) . In particular, in the approximation (|41|) . a fit to the 
CLEO and L3 data on the rj' — 7 transition form factor for Q 2 larger than 2 GeV 2 was 
recently undertaken in Ref. 0, yielding 

A 2 (l GeV 2 ) = -0.08 ±0.04, B 2 (l GeV 2 ) = 9 ± 12, (79) 

where, in the analysis in Ref. the initial scale in the evolution of the Gegenbauer 
moments was taken as /Xq = 1 GeV 2 . The estimates (179)1 can be translated in terms of 
the universal free parameters (Ho) and B^i/J^) (the Gegenbauer coefficients), fixed 
at the same initial scale /Xq of the DA evolution, yielding: 

Bi q \l GeV 2 ) = 0.02 ± 0.17, B ( 2 9 \l GeV 2 ) = 9.0 ± 11.5. (80) 
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We note that this analysis of the ?/ — 7 transition form factor yields an order of magnitude 
uncertainty in these parameters, taking into account the ±ler error. In addition, the 
Gegenbauer moments A 2 (l GeV 2 ) and -82(1 GeV 2 ) are strongly correlated and we refer 
to Fig. 3 of Ref. |3] where the correlation between these quantities is presented. 

The other process which allows to get an independent information on the Gegenbauer 
coefficients in the r]'g*g vertex is the inclusive decay T(IS') — > rj'X. Recently the rf- 
meson energy spectrum was measured in this process by the CLEO collaboration and 
presented in Ref. jH]. The CLEO data prefers the perturbative-QCD motivated form of 
the rj' — g transition form factor (^) for the hard part of the 7/-meson energy spectrum, 
Eff > 0.35Mt, where My = 9.46 GeV is the mass of the T(15')-resonance. Based on this 
observation, the form (|7Hj) for the function H(qf, 0, m 2 ,) resulting from the hard-scattering 
perturbative-QCD approach for the t]'g*g effective vertex function was adopted by us in 
Ref. jH] , obtaining independent constraints on the Gegenbauer coefficients. Unfortunately, 
the CLEO data on the decay T(IS') — > rj'X is statistically very uncertain for the end 
part of the //-meson energy spectrum, leaving large uncertainties in the determination of 
the Gegenbauer coefficients from this process alone. The combined fit to the data from 
the rf — 7 transition form factor and the T(15) — > rj'X decay leads to more restrictive 
constraints on these coefficients which we have presented and discussed in detail in Ref. [5] . 
The combined best fit of the Gegenbauer coefficients and the Gegenbauer moments yields 
the following (± lcr) values, respectively (SJ: 

B { 2 q) (til) = -0.008 ± 0.054, B { 2 9) (fi 2 Q ) = 4.6 ± 2.5, (81) 
A 2 (fj 2 ) = -0.054 ± 0.029, B 2 (fi 2 ) = 4.6 ± 2.7, 

where the starting point of the evolution in the //-meson DAs is taken as /z 2 , = 2 GeV 2 . 

The predictions for the rj'g*g effective vertex function F 71 > g * g (q 2 , 0, m 2 ,), with the second 
gluon on the mass shell (q 2 = 0), are presented in Fig. El for the time-like (left frame) and 
space-like (right frame) virtuality q\ of the off-shell gluon. The labels on the curves C, 
LL, and UR correspond to the central, lower-left, and upper-right points of the combined 
best fit of the Gegenbauer coefficients presented in Fig. 5 of Ref. jSJ; their values can 
also be read from Eq. (jHTj) . We note that the effect of including the ?/-meson mass in 
the vertex function is significant; it becomes crucial for the time-like gluon virtuality in 
the region q\ < 5 GeV 2 for the upper-right part of the combined best-fit region of the 
Gegenbauer coefficients (the curves labeled as UR), but also to a lesser extent for the 
central values of the fit parameters (the curves labeled as C). For the space-like gluon 
virtuality the 7/-meson mass effect is numerically not so strong. Nevertheless, it decreases 
the absolute value of the r]'g*g transition form factor by approximately 10 percent in the 
region of small gluon virtualities (\q 2 \ < 5 GeV 2 ) for the Gegenbauer coefficients from the 
lower-left part of the combined best-fit region (denoted by the curves labeled as LL) and 
for the central values of the fit (the curves labeled C). 

Summarizing this section, the results presented in Eqs. (J77J) and (J78J) (to be read with 
Eq. (JTJ)) are our principal analytic results for the vertex functions F v ig* g * (q 2 , q 2 , m 2 ,) and 
F ri /g* g ( y qf,0,m 2 ,), respectively. The vertex function F n ig* g (q 2 , 0, m 2 ,) is displayed numer- 
ically in Fig. |S] in the space-like and time-like regions of the gluon virtuality, based on 
perturbative QCD and our current knowledge of the Gegenbauer coefficients. Inclusion 
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Figure 5: The r]'g*g effective vertex function for the time-like (left frame) and space-like 
(right frame) gluon virtuality q\ when the second gluon is on the mass shell {q\ = 0). 
The solid and dashed curves are plotted for the vertex function with and without taking 
into account the ?/-meson mass, respectively. The labels C, LL, and UR correspond to 
the central, lower-left and upper-right points of the combined best fit of the Gegenbauer 



coefficients 



(/io) and B2 {fJ>o) presented in Fig. 5 of Ref. P and given in Eq. (JHTj) . 



of the 7/-meson mass effect has led to a rather reliable estimate of this vertex function 
(or the transition form factor) for q\>2 GeV 2 for the time-like region, with a power-like 
fall-off with qf, similar to the one seen in the electromagnetic transition form factors of 
the pseudoscalar mesons [2211211] • For the space-like gluon virtuality, the vertex function 
has a power-like fall-off as well, but the current uncertainty in the Gegenbauer coefficients 
still leaves a rather large dispersion in the vertex function, in particular for the region 
\Qi\ ~ 3 GeV 2 . We expect that more precise data, such as from the decay T(1S) — ► rj'X, 
will significantly reduce this parametric uncertainty. In the next section, we show that 
imposing the anomaly constraint on the vertex function F ri i g *g(qf,0,m^,) at q\ = con- 
siderably reduces this uncertainty for low values of q\ in the space-like region. 



5 An Interpolating Formula for the rfg*g Vertex Func- 
tion for Space-Like Gluon Virtuality 

As noted earlier, the QCD anomaly determines the value of the vertex F r} > gg (Q, 0, m 2 ,) 
for on-shell gluons. Denoting this by F£ , one has the following expression for this 
quantity [2]: 

F v % 9 = -4™ S K) • (82) 

For large off-shellness of the gluons, the form of the vertex function is determined by the 
perturbative QCD and is presented in Eq. (|77jl. We would like to write down an expres- 
sion for the vertex function in question which interpolates between the non-perturbative 
result ([E2J)> applicable at q\ = q\ = 0, and the perturbative QCD result (J77J), which holds 
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for large virtualities of the gluons. Such a formula is of considerable phenomenological 
interest. 

Let us concentrate on the case of the rj'g*g effective vertex function with one gluon (the 
second gluon, for defmiteness) being on the mass shell. Comparing Eqs. (Jfi7|) and (J7HJ) 
with the form (|T|. it is seen that both the quark- ant iquark and the gluonic parts of 
the rj' — g transition form factor have a pole singularity for the time-like virtuality at 
q\ = rat . However, there is no singularity for the space-like region of the gluon virtuality 
and the vertex function F v i g *g(qf,0,m 2 ,) is a smooth function of qf in this region (albeit 
numerically uncertain due to the imprecise knowledge of the Gegenbauer coefficients). 
It is also known from our earlier work [2] that the singularity for the time-like region 
can only be removed by including the transverse momentum effects in the 7/-meson wave- 
function, using the Sudakov resummation technique. Since we are ignoring the transverse- 
momentum effects in calculating the vertex function in question in this paper, we shall 
restrict ourselves to the interpolating function only in the space-like region, for which the 
transverse- momentum effects are known to be numerically small |2]. 

To that end, we work with the function H(qf, 0, rat) presented in Eq. (|75|l. It should 

be noted that the function G 2 f (1, rj) entering in H(qf, 0, rat) is very close to its asymptotic 
value 1/6 already at q\ ~ — 1 GeV 2 in the space-like region of the gluon virtuality (see 
Fig. HJ). So, to a very good approximation, the function G 2 , (1 , 77) can be replaced by its 
asymptotic value: 



47ra,(Q 2 



rat 



l + A 2 (Q 2 )-^B 2 (Q 2 ) 



(83) 



with the corresponding vertex function now given by F as (qf) = m^,H as (qf)/(qf — ra 2 ,). 
Thus, the dependence of H^q 2 ) on qf is coming only through Q 2 = \qf \ + ra 2 ,. Defined in 
this way, the function if as (qf ) is symmetric under the change q\ — > —qf, and has formally 
the following limit for on-shell gluons (qf = 0): 



Ana s ( rat) 

#as(0) = V^^/V 

rat 



l + A 2 (m 2 l ,)-—B 2 (m 2 i ; / 



34) 



Not unexpectedly, there is a substantial mismatch between the correct value of the vertex 
function for the on-shell gluons Ffi , as determined by the QCD anomaly, and the one 
obtained from the formal limit of the perturbative expression for the vertex function 
F v i g *g(0, 0, rat) = — if as (0), given in Eq. (|5^|i. To see this quantitatively, we study the 
numerical consistency of the two expressions, which yields the following condition on the 
Gegenbauer moments A 2 (rat) and B 2 (ral,): 



2VSn 2 f 2 



m 



A, 



ra r 



— B 2 (mi 
36 V " 



35) 



For the values of the Gegenbauer coefficients given in Eq. (JHTJ), this equality is badly 
violated. For example, for the set of the Gegenbauer moments A 2 (ra 2 )l ) = —0.11 and 
B 2 (ra 2 ,) = 2.87, corresponding to the coefficients called LL (yielding the largest value for 
the l.h.s. in Eq. (|85jl ). the l.h.s. in the above equation is about 0.43. For other allowed 
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values of the Gegenbauer coefficients, the mismatch is much more pronounced. This 
implies the presence of large non-perturbative contributions to the moments A 2 (m 2 ,) and 
However, it is certain that non-perturbative corrections in the vertex function 
F v 'g* g (qf, 0, ml,) are to be included over a larger region of qf. 

To model these non-perturbative effects, we propose the following modification of the 
perturbative result (JHHj) in the region of the small space-like virtualities q\ in terms of the 
function H(qf): 



H(qf) = H^qf) + [H A - ff M (0)] exp 



s 2 



(86) 



with the corresponding vertex function defined as: F(qf) = m 2 ; H (qf) / (qf — m 2 ,). Here, 
Ha = —F^gg is defined by the anomaly value for the r/'gg vertex in Eq. (|H2^1 and C s is an 
arbitrary dimensionless parameter. 

This formula gives, by construction, the correct value for the rj'gg vertex function 
with the on-shell gluons (i.e., for qf — 0), reproducing the QCD anomaly, and yields the 
correct perturbative-QCD behaviour for large values of qf. Depending on the value of 
the parameter C s , the perturbative result may set in rather fast. For small values of the 
gluon virtuality, expanding the exponential factor in qf/mf^, gives power corrections. In 
the numerical analysis presented below we take C s = 2, as it allows a smooth interpolation 
between the anomaly (at qf = 0) and the perturbative result for \qf\ > 2 GeV 2 . However, 
we must admit that this interpolating formula is by no means unique, and has to be 
checked against experimental data on processes involving the //-meson, or else compared 
with the results obtained using non-perturbative techniques, such as the lattice-QCD. 

The perturbative vertex function —F^qf) (dashed curves) and the interpolating func- 
tion —F(qf) (solid curves) are plotted as functions of qf in Fig. El (left frame). The function 
H as (qf) (dashed curves), defined in Eq. (JHHj) . and the corresponding interpolating function 
H(qf) (solid curves), introduced in Eq. ()86|) . are shown in the right frame of this figure. 
The curves are marked as C, LL, and UR, corresponding to the three sets of values for the 
Gegenbauer coefficients as in Fig. It is obvious that imposing the anomaly condition 
as a normalization point changes the t]'g*g vertex function in a significant way for low 
space-like values. For our choice of the interpolating function, this change is marked for 
\qf\ < 1 GeV 2 , increasing the absolute value of the function, and reducing the theoretical 
dispersion on the vertex function in this region arising due to the imprecise knowledge of 
the Gegenbauer coefficients in the perturbative expression. 



6 Conclusions 

We have calculated the rj'g*g^ effective vertex function in the perturbative QCD approach 
using the light-cone DAs for the r/'-meson with the inclusion of the 7/-meson mass. It is 
shown that if one of the gluons is on the mass shell, the pole-like behaviour of the 
i]' - gluon transition form factor emerges in this approach for both the quark-antiquark 
and the gluonic parts of the form factor, and the corresponding function H(qf,0, m^,) 
is perturbatively calculated. The Gegenbauer coefficients, required for a quantitative 
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Figure 6: The rj' — g transition form factor in the perturbative QCD approach and using 
the interpolating formula for the space-like region of the gluon virtuality q\. The left 
frame shows the functions -F a s(<?i) (dashed curves) and F{q\) (solid curves). The right 
frame shows the functions if as (g 2 ) (dashed curves) and H(qf) (solid curves). The labels 
on the curves are the same as in Fig. 03 



analysis of the vertex function, have been taken from the combined analysis of the 7/ — 7 
transition form factor and the T(IS') — > rj'X decay, reported by us earlier in Ref. 0. The 
corrections due to the 7/-meson mass are analyzed numerically, with the result that they 
are important for lower values of the gluon virtuality, in particular in the time-like region. 
An interpolating formula connecting the QCD anomaly value and the perturbative-QCD 
behaviour of the r] 1 - gluon transition form factor is presented for the space-like gluon 
virtuality, taking into account the ?/-meson mass, which modifies the vertex function 
significantly in the region \q\\ < 1 GeV 2 and reduces the theoretical dispersion in low \q\\ 
region considerably. 
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